Random Vortex-Street Model for a Self-Similar Plane Turbulent Jet 
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We ask what determines the (small) angle of turbulent jets. To answer this question we first 
construct a deterministic vortex-street model representing the large scale structure in a self-similar 
plane turbulent jet. Without adjustable parameters the model reproduces the mean velocity profiles 
and the transverse positions of the large scale structures, including their mean sweeping velocities, 
in a quantitative agreement with experiments. Nevertheless the exact self similar arrangement of 
the vortices (or any other deterministic model) necessarily leads to a collapse of the jet angle. 
The observed (small) angle results from a competition between vortex sweeping tending to strongly 
collapse the jet and randomness in the vortex structure, with the latter resulting in a weak spreading 
of the jet. 
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Introduction: We address the apparent universality of 
the self-similar structure of plain turbulent free jet. It 
had been amply documented[l~-7] that such jets contain 
dynamically dominant large-scale self-organized motions 
consisting of two lines of large-scale vortices centered at 
staggered position on the two sides of the jet: y{x) — 
±ab(x) , b{x) — ax, a ~ 0.85, where x is a stream- 
wise coordinate, b{x) is the jet half-width [8] and a is the 
jet angle. This structure carries about 75% of the total 
kinetic energy [7]. Importantly, the remaining turbulent 
energy randomizes in part the positions and amplitudes 
of the coherent vortices. 

In this Letter we first construct a deterministic 
vortex-street model that excellently reproduces the 
experimentally observed longitudinal mean velocity 
profile without any adjustable parameter. The model 
predicts the transverse positions of the large scale 
structures, including their mean sweeping velocities, in a 
quantitative agreement with experiments. Nevertheless 
both the existence of a regular array of vortices and the 
turbulent randomization are crucial for understanding 
the (relatively small) observed angle a, 



cxp 
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the jet spreading. Therefore, secondly we develop a 
random-vortex street model, that leads to an estimate 
of the jet angle amod ~ 0.13 (reasonably close to the 
observed angle) as a result of a competition between the 
regular and the stochastic components. 

The deterministic model: To model the coherent 
structure of a plane jet we employ vortices with a finite 
core size [3] . Here we choose a simple algebraic model for 
the vorticity uj{r) and angular velocity Vip{r) that allows 
an analytic treatment: 
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where r is the distance from the vortex center. To charac- 
terize the self-similar structure of the jet we enumerate 



the vortices by their sequential number — oo < n < oo 
and introduce a scaling parameter A > 1 such that the 
x-position of vortex centers are Xn — a;oA", where xq 
is the stream-wise position of a reference vortex. The 
y-position of the vortices is chosen at the points of inflec- 
tion of the mean velocity profile where the production 
of vorticity has its maximum: ?/„ = (— l)"aQ;a;„, where 
both a and a need to be determined below. 

Note that we construct the jet exactly self-similar; in 
experiments this is achieved only asymptotically away 
from the orifice, about 30 times the width of the orifice. 
To flush out the scaling properties of the mean stream- 
wise velocity Vx{x, y) and the velocity circulation r„ no- 
tice that for the free jet the stream-wise momentum flux 
M is a;- independent: 
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Since b{x) scales as x, for constant M the centerline ve- 
locity Vx{x,0) should necessarily scale like l/y^. Simi- 
larly, one estimates r„ ~ V^(x„, 0)6(a;„) oc y/x^ oc \/A". 
Therefore r„ = -f{~V\y\ with some dimensional con- 
stant 7. 

Note that this model, if used as initial condition for 
the Euler equation, will lose its self-similarity immedi- 
ately. Dynamically suitable model that reflects the self 
similarity of the observed jet will call for a very very 
much more detailed and parameterized description. The 
analytic simplicity of the present model results from the 
fact that the self-similarity should be interpreted statis- 
tically rather than dynamically; in the real jet situation 
we must allow for the coalescence or sometime creation 
or death of vortices. 

At an arbitrary position x, y and time t the stream- 



wise velocity field, produced by all the vortices, is: 

U.ix.y) = -L^_(_VArr„/i?2, (3a) 
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where the core radius of the n-th vortex c„ = CqA". Sim- 
ilarly one gets an equation for the transverse velocity 
Uy{x,y) by replacing —Yn Xn in the numerator of 
Eq. (3). 

We define the mean (in time) velocity profiles, 
Vx{x,y) = {Ux{x,y)) and Vy{x,y) = {Uy{x,y)) and the 
normalized coordinates ^ = x/b{x)X = y/h{x). As the 
next step we define a normalized self-similar mean veloc- 
ity, which is x-independent in the self-similar regime: 



V^{(:) = V{x,y/b{x))/V,{x,Q) 
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Since the actual experimental value of the spreading 
angle ctoxp is small, and since we will find below that the 
profile of the longitudinal mean velocity is very weakly 
dependent on a, it is very worthwhile to proceed and an- 
alyze precisely the limit a — *■ 0, where the vortices are 
arranged on two parallel lines y = ±a, separated in the 
streamvise direction by d. There we can simplify dramat- 
ically the calculation of the mean stream-wise velocity 
from Eq. (3), in the form of a single integral. For a — > 
we write A" 1 + (A — l)n and substitute in Eq. (3) 
Xn = xq + ndh, d ^ {\ - l)xo, yn = (-l)"a6, c„ = cb, 
where b is the x-independent jet half-width: 
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The average with respect to time was replaced by aver- 
aging in the ^ direction over one period of oscillations 
(between two consecutive vortices at — c? and d). In each 
term in the sum (separately for odd and even n) we can 
change the integration variables ^' = S_ + 2dn and even- 
tually collapse the whole sum into two integrals, finding 
ioYV.,{bQ = {U^[x,y)): 
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Notice that Vy{C,) vanishes in the hmit a ^ due to sym- 
metry. While these resuhs were developed for the limit 
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FIG. 1: Color online. Profiles of the normalized stream-wise 
mean velocity: Black squares - experimental points [7], black 
dashed line - profile (9) obtained from the Prandtl closure 
with constant turbulent viscosity approximation, red dash- 
dotted line line - Eq. (7) (a ^ limit) with a = 0.747, 
c « 1.1. Black solid line - profile (8) with a = 0.1 and the 
same values of a and c. 



a — !■ 0, we will demonstrate that they pertain excellently 
well also for small a. The normalized profile according 
to Eqs. (4) and (6a): 

v|(C)^[/+(C) + /-(C)]/[/+(o) + /-(o)] . (7) 

This profile depends only on a and c. With the definition 
of the width of the jet, b{x), we now demand Vj;(l) = 
^, giving us one relation between a and c. The second 
relation between these parameters is determined by the 
position of the inflection point: (i^V|(C)/dC^ = at C = 
a. Solving the two conditions together we find a ~ 0.747 
and c = 1.1. The resulting profile (7) with these values 
of the parameters is shown in Fig. 1. This profile V|(C) 
obtained for a = should be compared to the profile 
V^{C) computed for finite a [9]: 



[j+(C) + J-(C)]/[J+(o) + J-(o)], 
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y/[{ax ± 0' + (X - 1)V«']' + C^X' 



In Fig. 1 we compare these two profiles and find that 
as expected they are practically indistinguishable for the 
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0.1. Moreover, the profile 

< 



experimental value a, 

Eq. (8) is almost insensitive to a for small a ^ 0.15. 
In the same figure we compare the model profile to the 
experimental data [7], and find excellent agreement. In 
contrast, the profile 



KL..(C)-cosh-2 ln(l + V2) 
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computed using the "uniform turbulent viscosity" as- 
sumption [10], deviates strongly from the data for ^ > 1. 

Moreover, the present model fits the data in more than 
one way. One can actually measure the transverse posi- 
tions yn of the large scale vortices in units of the 6(x„). 
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FIG. 2: Vortex structure (a fragment from experimental 
proper orthogonal decomposition, Ref. [7]). 



relative success of the deterministic vortex-street model 
was possible because fluctuations of the vortex param- 
eters have little consequences on the mean values; they 
can be neglected in zeroth-order approximation in de- 
scribing the mean velocity profiles. This is not the case 
for the spreading angle a, which tends to collapse in the 
deterministic limit. 

To find the spreading angle a we should consider the 
normalized velocity fluctuations in the transverse direc- 
tion y, 



In the self-similar regime of the jet the ratio yn/b{xn) is 
n-independent and best estimates [7] give approximately 
0.85 in a reasonable agreement with our value of a « 0.75. 
Finally, the measured value of the stream-wise velocity 
of the coherent structures in units of the centerline ve- 
locity is about 0.60 — 0.65 [4, 5]. Our model predicts 
V|(a) « 0.57. We remind the reader that all this agree- 
ment is achieved without using any experimental data. 

Having modeled so successfully the stream- wise veloc- 
ity, we can now evaluate the transverse velocity by using 
the incompressibility constraint and Eq. (4). We derive 
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For small a one can neglect the a dependence of the 
integrand in the right hand side of Eq. (10), replacing Vj 
by its Of ^ limit V|. Now we can compute the ratio 
a' = Vy{x,y)/Vx{x,y) dX y — ab{x), finding that the 
sweeping angle of the vortices a' « 0.3 a is considerably 
smaller than the angle aa « 0.75a that is required to 
preserve the angle self-similarly. This result means that 
the mean velocity of the vortices is not maintaining the 
self-similar structure at any (small) value of a, tending 
to collapse the structure to two parallel lines. We expect 
this result to be generic for any deterministic model. 
Random Vortex-Street Model. Now we have to ac- 
count for the fact that the real jet does not display ideal 
vortex-street structure, see, e.g. Fig. 2. The parame- 
ters of individual vortices fluctuate, and we also see the 
coalescence of vortices as mentioned above. Moreover, 
pairs of vortices have finite life-time, as we see in Fig. 2, 
that two vortices are about to disappear at Uct « —8 6. 
We reiterate that the finiteness of the vortex life-time is 
required by the self- similarity, i.e. the sweeping velocity 
is oriented along the line that connects the vortex cen- 
ters. Indeed, the sweeping velocity scales like l/\/x, the 
vortex density scales like 1/x, therefore the flux of the 
vortex number scales like 1/x^^^ i.e. decreases with x. 
This means that vortices disappear during their sweeping 
downstream, i.e. they have finite life time, due to their 
coalescence, or decay, or whatever. It means, as said 
above, that self-similar jet structure must be understood 
in the statistical sense and randomness of the vortex pa- 
rameters is a necessarily element of the jet structure. The 
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One simple way to introduce randomness in the model 
is to account for the fluctuations of the core size of the 
vortices, c„ = c {l + 5n) with uncorrelated statistical fluc- 
tuations: (S„ = , dn^n' = 5'^^n,n'- In the limit a — > 
we replace the parameter c in 7?,^ by c„ and (for small 
fluctuations, 5^ ^ 1 ) expand Eqs. (5) with respect of 
5n ■ Averaging the result with respect to the fluctuations 
and over longitudinal position in the jet we get: 
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In accordance with these results we present K as a sum 
of the regular part K'^'^^, that originates from x-periodic 
velocity fluctuations, and the random part K'^^'^, that 
originates in our model from statistical fluctuations of the 



core radius: K 
and K''''' E 
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related in Eqs. (12) with E^^ |, and E^^'^,. Experimental 
information about K is provided in Fig. 3. We see from 
the data that on the axis of the jet « 0.04 and at the 
centers of the vortices K^{( — a) ^ 0.03. We use this 
data to fix the parameters of the model with the results 
d = 1.45 and = 0.03. With these parameters the 
regular contribution K'^°s(a) w 0.024 while the random 
contribution is considerably smaller, K^^^{a) « 0.006. 

In this way we have separated the regular contribution 
to K that does not effect the spreading angle, from the 
random one, K'^^'^, that leads to the jet spreading by 
turbulent diffusion with a characteristic velocity y/K™^. 
We thus estimate 
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FIG. 3: Experimental profile [7] of the mean square of the 
cross-stream component of the fluctuating velocity (dots con- 
nected by line). In a continuous line we present the profile 
predicted by the random model. The dashed line is the profile 
of the regular contribution. The difference at a is the source 
of spreading of the jet, cf. Eq. (13). 

This estimate of the tangent of the angle of the jet 
is slightly higher than the experimental value (about 
0.1). We note however that we made no effort to model 
quantitatively the random structure of the vortices. For 
example, looking at Fig. 2 we note that the vortices are 
elliptic rather than round, this will definitely contribute 
to lowering the random part of K^. One could input 
this knowledge, with a price of additional fit parameters. 
We submit to the reader that this is not the goal of 
the present calculation, we are not interested here in 
a quantitative model of a jet spreading. We wanted 
to understand what is the physical reason for this 
phenomenon and why the jet angle is small. We believe 
that the model worked out above is sufficient; our answer 
is that any deterministic model will necessarily close the 
jet angle, allowing only a self similar parallel jet. The 
reason for the opening angle is randomness, which is 
the flip side of turbulence which cannot be avoided even 
when 75% of the energy is in the coherent structure. The 
randomness contributes partially to the value , which 
is rather small per se. Therefore (without which 

the opc;ning angle will tend to zero) is even smaller and 
thus resulting spreading angle (13) is indeed small. 

Summary. In this Letter we offered a simple model of a 
plane jet with the aim of understanding the fundamental 
existence of a statistically self similar structure with a 
small opening angle. The deterministic part of the model 
excellently reproduces the experimental profiles of the 
mean velocity without using any adjustable parameter. 
For an opening angle we must introduce randomness. In 



the context of the present simple model we ascribed the 
randomness to the core size, and demonstrated that this 
is sufficient to provide an opening angle of the correct 
order of magnitude. Needless to say we do not pretend 
that this model describes correctly the full randomness 
which is due to turbulence, but only underlines the 
crucial role of the random components in opening up the 
jet. One can provide a more precise model with better 
agreement with the measured angle and the profiles of 
second order quantities, but this calls for additional 
adjustable parameters. Such a detailed model is not the 
aim of this Letter which concentrated on understanding 
the basic physics of the phenomenon. 
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